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0\ 1 Abstract 
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Using a single trial function we perform an accurate calculation of the ground 
state \o g of the hydrogenic molecular ion in a constant uniform magnetic 
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field ranging — 10 13 G. We show that this trial function also makes it possible 
to study the negative parity ground state la u . We obtain that over the whole 
CN ■ range of magnetic fields studied, the calculated binding energies are in most 

q ■ cases larger than binding energies obtained previously by other authors using 

I/"") ■ different methods. 
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I. INTRODUCTION 



For a long time the behavior of atomic and molecular systems in a strong magnetic field 
of order B > 10 9 G has attracted considerable attention. The main interest comes from both 
astrophysics - studies of white dwarfs and neutron stars (see 0,@], and also, for example 
0HJ], and the review [||] and references therein), as well as from chemistry - formation of 
unusual chemical compounds whose existence is impossible without strong magnetic fields 
H (for a review, see, for example, J7J and references therein). 

There are many studies of the hydrogen atom H - the simplest atomic system - in a 
strong magnetic field, while the hydrogen molecular ion H£ - the simplest molecular system 
which is stable with respect to dissociation H% — > H + p - is much less explored. One of the 
major drawbacks of many of these studies is a restricted domain of applicability: they are 
accurate in the weak magnetic field region but are inappropriate for the strong magnetic field 
region and vice versa. The goal of the present Note is to carry out an accurate variational 
calculation of H 2 in magnetic fields ranging from up to 10 13 G [] using a unique simple 
trial function equally applicable for any value of the magnetic field strength. We restrict 
our consideration to the case where the magnetic field is directed along the axis of the 
molecule, which is evidently the optimal configuration leading to the lowest energy. Our 
main perception is that the calculations should not be technically complicated and also 
easily reproduced, while the trial function should be simple enough to allow further analytic 
and numerical investigations. 



II. CHOICE OF TRIAL FUNCTIONS 

A constructive criterion for an adequate choice of trial function was formulated in || 
and further development was presented in P,|T0||. In the simplest form the criterion is the 
following. The trial function ^t{x) should contain all symmetry properties of the problem 
in hand. If the ground state is studied, the trial function should not vanish inside of the 
domain where the problem is defined. The potential V t (x) = 9 , for which the trial 
function is an eigenfunction, should reproduce the original potential near singularities and 
also its asymptotic behavior. The use of this simplest possible recipe has led to a unique 
one-parameter trial function, which, in particular made it possible to carry out the first 
qualitative study of the ground state of the hydrogen molecule H 2 in the region of both 
weak and strong magnetic fields ||11|| . Later a few-parameter trial function was proposed 



for a description of the hydrogen atom in an arbitrary magnetic field, which led, for the 



low-excited states, to an accuracy comparable with the best calculations fll0| , |12 l. 

Now we wish to apply the above-mentioned recipe to the ion H 2 ■ We work in the Born- 
Oppenheimer approximation. Let us first introduce notation (see Fig.l). We consider two 
attractive centers of charge Z situated on the z-axis at a distance R/2 symmetrically with 
respect to the origin. The magnetic field of the strength B is directed along the z axis and 



x where the relativistic corrections can still be neglected (see a discussion in || and references 
therein) 
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r\ t 2 are the distances from the electron to the first (second) center, respectively. The quantity 
p is the distance from the electron to the z-axis. Through the paper the Rydberg is used 
as the energy unit. For the other quantities standard atomic units are used. The potential 
corresponding to the problem we study is given by 

2£ ! _2Z_2Z BV (1) 
R r, r 2 4 ' 1 ' 

where the first term has the meaning of the classical Coulomb energy of the interaction of 
two charged centers. 

One of the simplest functions satisfying the above recipe is the Heitler-London function 
multiplied by the lowest Landau orbital: 

q, 1 — e -aiZ(r 1 +r 2 )-/3 1 Bp 2 /4^ ^) 

where a±, Pi are variational parameters. It has in total 3 variational parameters if we include 
the internuclear distance R, in the search for the equilibrium distance. It is known that 
in the absence of a magnetic field a function of the Heitler-London type gives an adequate 
description of diatomic systems near their equilibrium position. The potential corresponding 
to this function is: 



V x = 2Z 2 a\ - Bfc - 2Za x (— + - 

Wi TV 



f3 2 B 2 p 2 



2 

+ 2Z 2 —^- 

r x r 2 



p 2 + (z- R/2) (z + R/2)] + Z ai Bp lP 2 (- + -)• (3) 



It is clear that this potential reproduces the original potential (1) near Coulomb singu- 
larities and at large distances, \x\ — > oo. 

The Hund-Mulliken function multiplied by the lowest Landau orbital is another possible 
trial function: 

^ 2 = ( e -^Z ri + e - aa Z n \ e -fhB?/i ^ ( 4 ) 



where 0:2, 02 are variational parameters. It is well known that this function, in the absence 
of a magnetic field, describes the region of large internuclear distances. The calculations we 
performed show that this property remains valid for all magnetic fields up to 10 13 G. Like 
(2), the trial function (4) is characterized by 3 variational parameters. In order to take into 
account both equilibrium and large distances, we should use an interpolation of (2) and (4). 
There are two natural ways to interpolate: 

(i) a non-linear superposition: 

^ _ ^ e -a3-^»*i-a4^r2 _|_ - e -«4^ri-a 3 Zr2~j e -/33-Bp 2 /4 ^ 

where 0:3,0:4, (3 3 are variational parameters. The parameter a = ±1 depends on 
whether positive parity la g or negative parity la u states we consider. The func- 
tion (5) is a modification of the Guillemin-Zener function used for the description of 
the molecular ion H^. If 0:3 = 0:4, the function (5) reduces to (2) and if o 3 = 0, it 
coincides with (4). In total there are 4 variational parameters characterizing the trial 
function (5); 
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(ii) a linear superposition of (2), (4) 



(6) 



where the relative weight of (2) and (4) in (6) is taken as an extra variational parameter. It 
is a 7-parameter trial function. 

Of course, as a natural continuation of the above interpolation procedure one can take 
a linear superposition of all three functions: the modified Heitler-London, Hund-Mulliken 
and Guillemin-Zener functions (2), (4), (5) 

#4 = 4j_l*3-l + ^3-2^3-2 = ^3-1^3-1 + ^1*1 + ^2*2 , (7) 

where again, as in the case of the function (6) the relative weights of different components 
are variational parameters. In total, the trial function (7) is characterized by 10 variational 
parameters. Most of our calculations will be carried out using this function. The minimiza- 
tion procedure was done using the standard minimization package MINUIT from CERN-LIB 
on a Pentium-Pro PC. 



III. RESULTS 

In order to present our results we begin with the study of the dependence of our vari- 
ational results on different trial functions (see Table 1). It turned out that all variational 
parameters are of the order of one independently on the value of the magnetic field strength. 
In the Table 2 we give a comparison of our calculations with the best known results. Since 
we are doing a variational study of the problem, the (binding) energies obtained represent 
(lower) upper bounds to the exact energies. No need to mention that many calculations 
were performed for the case of the H£ ion in absence of a magnetic field. Our results are in 
agreement with the best calculations within an absolute accuracy of 10~ 5 . For all studied 
values of the magnetic field (< 10 13 G) our results for binding energies exceed the best known 
results at the present. Among all previously made calculations, we should emphasize that 



those performed by Wille |]T3[ have the most extended domain of applicability. The accuracy 
of these results is almost as good as the accuracy of our results for magnetic fields B < 10 11 
G. However for bigger magnetic fields the accuracy of his results fells down drastically. 

Fig. 2 shows the electronic density distribution as a function of magnetic field. For a small 
magnetic field the distribution has two clear maxima corresponding to the positions of the 
centers. The situation changes drastically for magnetic fields of the order of B ~ 3 x 10 11 
G, where the probability of finding the electron in any position between the two centers is 
practically the same. For larger magnetic fields the electron is preferably located between 
the two centers with the maximum in the middle, z = (see Fig 2 (e), (f)). Due to a loss of 
accuracy this phenomenon was not observed in |L4| . It is worth noticing that for all magnetic 



fields studied the region of large internuclear distances is dominated by the Hund-Mulliken 
function (4). 

In Table 3 the results of the study of the ground state of negative parity (la u ) are 
presented. Let us note that in the absence of a magnetic field the electron is not bounded 



possessing a shallow minimum (see, for example, []T3| . p~5|] ) . However, in a magnetic field the 



electron becomes bounded in agreement with general expectations || (see also [jig] ). 
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In conclusion we want to emphasize that we present the most accurate calculations for 
the ground state energies and equilibrium distances of the molecular ion in magnetic 
field. Unlike the majority of other studies our calculations stem from a unique framework 
covering both weak and strong magnetic field regimes. 
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(a) 



(b) 




FIG. 2. The electronic probability density of the ground state la g at y = as a function of the 
magnetic field strength, (a) B = 0, (b) B = 10 9 G, (c) B = 10 10 G, (d) B = 10 11 G. 
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TABLES 



TABLE I. Ground state energies and the equilibrium distances for the ground state la g of the 
ion H2 obtained using different trial wave functions as a function of the magnetic field strength. 

B ^1 # 2 *3-l ^3-2 *4 



E 

R 



-1.16277 
1.84678 



-1.17299 
2.00349 



-1.20488 
1.99799 



-1.20488 
1.99810 



-1.20525 
1.99706 



10 9 G 



E 
R 



-1.11760 
1.80967 



-1.11849 
1.93353 



-1.15007 
1.92378 



-1.15013 
1.92375 



-1.15071 
1.92333 



10 10 G 



E 
R 



1.11835 
1.21434 



1.15713 
1.24328 



1.10098 
1.25042 



1.10095 
1.25013 



1.08993 
1.24640 



10 11 G 



E 
R 



35.1065 
0.59575 



35.2351 
0.59206 



35.1016 
0.60065 



35.0750 
0.58954 



35.0374 
0.59203 



10 12 G 



E 
R 



408.539 
0.28928 



408.837 
0.28372 



408.539 
0.28927 



408.424 
0.28422 



408.300 
0.28333 
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TABLE II. Comparison of present calculations for the ground state la g of equilibrium distance 
R, energy E and binding energy BE = B — E with the results of other calculations [4,13-17]. 

B (Gauss) R (a.u.) E (a.u.) BE (a.u.) Source 



B = 



B = 10 9 G 



B = 10 10 G 



B = 10 11 G 



1.9971 
2.0000 
1.997 
1.997 
2.0000 

1.9233 
1.924 
1.921 
1.90 

1.2464 
1.246 
1.159 
1.10 

0.593 
0.593 
0.62 



-1.20525 
-1.20527 
-1.20527 
-1.20526 
-1.205268 

1.15072 
1.15072 



1.08979 
1.09031 



35.0362 
35.0428 



1.57616 
1.57616 
1.5757 
1.5529 

3.1646 
3.1641 
3.0036 
3.0411 

7.5080 
7.5013 
7.35 



Present 
Teller 



Wille |L4| 
Peek-Katriel |]| 
Wind pj 



Present 



Wille g4| 
Peek-Katriel H 
Lai-Suen []l7| 

Present 



Wille g4| 
Peek-Katriel @ 
Lai-Suen UTTf] 

Present 
Wille |14| 
Lai-Salpeter M 



B = 5 x 10 11 G 



0.351 
0.350 
0.35 



199.238 
199.315 



13.483 
13.406 
13.38 



Present 



Wille |4| 
Lai-Salpeter 



B = 10 12 G 



0.283 
0.278 
0.28 



408.300 
408.566 



17.141 
16.875 
17.06 



Present 
Wille g. 
Lai-Salpeter 



B = 2 x 10 12 G 



0.230 
0.23 



829.274 



21.609 
21.54 



Present 
Lai-Salpeter 



B = 5x 10 12 G 



0.177 
0.18 



2098.3 



28.954 
28.90 



Present 
Lai-Salpeter 



B = 10 13 G 



0.147 
0.15 



4218.7 



35.752 
35.74 



Present 
Lai-Salpeter 
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TABLE III. Energy, binding energy and equilibrium distance for the ground state of of 
negative parity la u . 



B (Gauss) 



R (a.u.) 



E (a.u.) 



BE (a.u.) 



Source 



B = 

B = 10 9 G 

B = 10 10 G 

B = 10 11 G 
B = 10 12 G 
B = 10 13 G 



12.746 
12.55 

11.039 
10.55 

6.4035 
4.18 

3.7391 

2.4329 

1.7532 



-1.00010 
-1.00012 

-0.92063 
-0.917134 

1.6585 
2.1294 

36.945 

413.92 

4232.6 



1.00010 
1.00012 

1.34608 



2.59592 

5.59901 
11.519 
21.851 



Present 

Peek-Katriel 111 



Present 

Peek-Katriel |1| 



Present 

Peek-Katriel pfl 



Present 
Present 
Present 
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